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Abstract 

We consider a class of partial differential equations with Carlitz derivatives over a 
local field of positive characteristic, for which an analog of the Cauchy problem is well- 
posed. Equations of such type correspond to quasi-holonomic modules over the ring of 
differential operators with Carlitz derivatives. The above class of equations includes some 
equations of hypergeometric type. Building on the work of Thakur, we develop his notion 
of the hypergeometric function of the first kind (whose parameters belonged initially to 
Z) in such a way that it becomes fully an object of the function field arithmetic, with the 
variable, parameters and values from the field of positive characteristic. 
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1 INTRODUCTION 

oo 

Let K be the field of formal Laurent series t = ^ S,jX^ with coefficients from the Galois 

j=N 

field Fg, .^AT 7^ if t 7^ 0, g = p" , v G Z_|., where p is a prime number. It is well known that 
any non-discrete locally compact field of characteristic p is isomorphic to such K. The absolute 
value on K is given by \t\ = , |0| = 0. This absolute value can be extended in a unique way 
onto the completion of an algebraic closure of K. See [H] for a detailed description of the 
extension procedure (valid for any non- Archimedean valued field). 
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An important feature of analysis over K and initiated by Carlitz ^ and developed 
by Wagner, Goss, Tliakur, the author, and many others (see surveys in [71 [21], [I6j) is the 
availability of many non-trivial Fg-linear functions, that is such functions / defined on Fg- 
subspaces Kq C with values from Kc, that 

/(tl+t2) = /(tl) + /(t2), /(«t)=«/(t), 

for any t,ti,t2 G Kq, a G ¥q. Such are, for example, polynomials and power series of the form 
(^kf^ ■ Within this class, there are analogs of the exponential and logarithm, the Bessel and 

hypergeometric functions, the polylogarithms and zeta function. 

The basic ingredients of the calculus of F^-linear functions [9] are the Frobenius operator 

Tu = u'^, the difference operator 

Au{t) = u{xt) - xu{t) (1.1) 

introduced in [3], and the nonlinear (Fg-linear) operator d = r^^A called the Carlitz derivative. 
The latter appears, as an analog of the classical derivative, in the theory of ordinary differential 
equations for Fg-linear functions, which has been developed both in the traditional analytic 
direction (the Cauchy problem [10, [13], regular singularity [TT], new special functions defined 
via differential equations [11]) and within various algebraic frameworks (analogs of the canonical 
commutation relations of quantum mechanics [8l,'9], umbral calculus [12], an analog of the Weyl 
algebra [TOl [T5l [2J). Note that in our situation the meaning of a polynomial coefficient in a 
differential equation is not a usual multiplication by a polynomial, but the action of a polynomial 
in the operator r (similarly, for holomorphic coefficients). Thus, the operator A = rd can be 

seen as the counterpart of t—. 

As an example, consider Thakur's hypergeometric function [191 1201 [21] (in this paper we 
deal only with what Thakur calls "the first analog of the hypergeometric function for a finite 
place of ¥q{xy). For n E Z+, a e Z, denote = [n][n - 1]" ■ ■ ■ [1]""'', [n] = x"" - x, 
L„ = [n][n -!]■■■ [1] (n > 1), Do = Lq = 1, 

Dn+a-1 , if a > 1; 

(«)n = < f-l)"-"L:f_„, if « < 0, n < -a; (1.2) 

if a < 0, 72 > —a. 

For ai,(3i G Z, such that the series below makes sense, Thakur's hypergeometric function is 
defined as 

Thakur has carried out a thorough investigation of the function (1.3) and obtained analogs 
of many classical properties. In particular, he found a differential equation for the function 
(1.3) (see [IH]). Using the commutation relation dr — rd = [1]^^'^/, where I is the identity 
operator, we can write the equation from [19] in the form 

n(A - [-a,]) - (^{A - [-/?,]) j rfj = (1.4) 
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convenient for us in a sequel. For r = 2, s = 1, this can be written in a familiar- looking form: 
the function u = 2-Fi(«, P', 1', z), a, P,'y E Z, satisfies the equation 



r(l - T)d\ + {([-1]^ + [-a] + [-f3]) T - [-7]} du - [-a][-/3]u = 0. 



(1.5) 



It must be remembered that everywhere in this theory we have to deal with nonlinear equations: 
the operators r and d are only Fg-linear, not K-, nor i^'c-hnear. 

A natural next step in developing analysis over K is to try to consider partial differential 
equations with Carlitz derivatives. However, here we encounter a serious difficulty noticed in 
|15] : the Carlitz derivatives with respect to different variables do not commute. Considering 
the Carlitz rings of "differential operators" in the above sense in P2|, the author found a 
class of partial differential operators (acting on an appropriate class of functions of several 
variables), which nevertheless possess reasonable properties. Operators from this class contain 
the derivative d in only one distinguished variable, and the linear operator A in every other 
variable. Such operators d and A (in different variables) do not commute either but satisfy a 
simple commutation relation. 

In this paper we pursue this idea further, introducing a class of partial differential equations 
of the above type, for which an analog of the Cauchy problem (the initial value problem 
for a differential equation) with respect to the distinguished variable is well-posed. In this 
regard, such equations may be seen as function field analogs of classical evolution equations 
of mathematical physics. On the other hand, it is easy to notice an analogy between our 
new equations and general hypergeometric equations of classical analysis (see, for example, 
[H [5] ; note also a different non- Archimedean generalization of the hypergeometric function [6] 
where the functions, defined on a non- Archimedean field, are complex- valued) . In this wider 
framework, it becomes clear how to extend the definition (1.3) of Thakur's hypergeometric 
function to the case of the parameters from Kc, so that the resulting function is locally analytic 
with respect to the parameters and coincides with Thakur's one (up to a change of variable), 
as the parameters are of the form [— n], n G N. Thus, now Thakur's hypergeometric function 
becomes fully an object of the function field arithmetic, with the variable, parameters, and 
values belonging to Kc- 

The structure of the paper is as follows. In Section 2, we prove the well-posedness of the 
Cauchy problem for the above class of partial differential equations. We show (Section 3) that 
these equations generate quasi- holonomic modules (in the sense of [IS]). Section 4 is devoted 
to the general hypergeometric function. Analogs of the contiguous relations for the latter are 
given in Section 5. 

2 Cauchy Problem 

Denote by J-'n+i {n > 1) the set of all germs of functions of the form 



where Cm,h,...,i„ G Kc are such that all the series are convergent on some neighbourhoods of the 
origin. Let J-'n+i be the set of polynomials from J-'n+i- 




m 



(2.1) 
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Below d will denote the Carlitz derivative in the variable z, while Aj will mean the difference 
operator (1.1) in the variable sj. In the action of each operator d,Aj on a function from J^n+i 
(acting in a single variable) other variables are treated as scalars. Obviously, linear operators 
Aj commute with multiplications by scalars: A^A = AA^, A G Kc, while d\ = X^^'^d. We have 
also the following commutation relations: 

dT-Td= [l]i/«J, dAj - Ajd = [lY/^d, AjT - tAj = [1]t (j = 1, . . . , n). (2.2) 

The Carlitz ring Sl^+i (see [I5j) is generated by r, d, Aj (j = 1, . . . , n) and scalars A G Kc, and 
the action of any operator from St^+i on J-'n+i is well-defined. 
Let us consider equations of the form 

{P(Ai, ...,An) + Q(Ai, . . . , An)d} u = (2.3) 

where P, Q are non-zero polynomials with coefficients from Kc- We look for a solution u G J^n+i 
of the form (2.1) satisfying an "initial condition" 

lim z~^u{z, Si, . . . , Sn) = Uo{Si, ...,Sn) (2.4) 

2— >0 

where mo(si, . . . , s„) is an F^-linear holomorphic function on a neighbourhood of the origin. 
The condition (2.4) (similar to the initial conditions for "ordinary" differential equations with 
Carlitz derivatives [10]) means actually that the coefficients co,ii,...,i„ of the solution (2.1) are 
prescribed for any ii, . . . , 

Below we use the notation [oo] = —x. Then [n] — > [oo], as n — > oo. 

Theorem 1. Suppose that 

Q{[ii], . . . ,[in]) 7^ for allii, ... ,in = 0,1, ... ,oo. (2.5) 
Then the Cauchy problem (2.3)-(2.4) has a unique solution u G J^n+i- 

Proof. It is easy to see that 

n ) |0, ifij = 0. 

The identity Dm = [m]D'^_-^ implies the relation 



di — ] = ( iim^O; 
Dm) [0, ifm = 0. 



Therefore for a function (2.1) we get 

CO oo min(ii,...,i„) 



*i =1 in=l m=l 



1/q gn-l 

m,ji,...,i,i 1 ■ ■ ■ n 



m—l 
oo min(ji,...,jr„) 



• • • ''i/+l,ii+l,...,j„+l'^l • • • ■^n 2^ ■ 

jl=0 i„=0 1^=0 " 



4 



Next, 



oo oo mm(ii,...,i„) 



il=0 in=0 m=0 



D. 



m 



Writing a similar formula for Q{Ai, . . . , A„)m and substituting all this into (2.3) we find 
that 

oo oo min(ji,...,i„) , 

41=0 «„=0 m=0 

€---€^-^ 

for arbitrary values of the variables. Hence, we come to the recursion 



q 

Cm+l,ii+l,...,in+l ~ ~^m,ii,...,ir. 



Since all the elements Co,ii,...,i„ (ii, . . . , = 0, 1, 2, . . .) are given, from (2.6) we find all the 
coefficients of (2.1). 

The set {[i], i = 0, 1, . . . , oo} is compact in K. Therefore the condition (2.5) imphes the 
inequality 

|Q([^i],...,W)| >/i>0 (2.7) 
for all ii, . . . , i„ = 0, 1, 2, . . . , oo. Note also that | [i] | = for any i, and 

|co,n,...,.J < Cr^"+-^^^" (2.8) 
for some positive constants C and r, since the series for the initial condition 



oo oo 



Uq{Si, . . . , S„) - ^ . . . ^ Co,ii,...,i„Si . . . 



n=0 in=0 

converges near the origin. 
By (2.7) and (2.8), 

|Ci,n+i,...,i„+i| ^ o^r 
(where Ci > does not depend on ii, . . . , i„), thus 

|C2,n+2,..,i„+2| S (^1 r , 

and, by induction 

|Q,il+i,...,j„+«| ^ "-^l 

for any / > 0. This means that for any ji, . . . ,jn > i, 

|Qjl,...Jn| S , O2 > U, 
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which, together with the equahty 

\Drn\ = q , 

imphes the convergence of the series in (2.1) near the origin. ■ 

Remark. It is easy to generahze Theorem 1 to the case of systems of equations, where P and 
Q are matrices whose elements are polynomials of Ai, . . . , A„. In this case, instead of (2.5) we 
have to require the invertibility of Q{[ii], . . . , [in]) for alHi, . . . , i„ = 0, 1, . . . , oo. In an obvious 
way, this generalization covers also the case of a scalar equation of a higher order in d. 

Below we will consider in detail some specific examples of the equation (2.3). However, 
before that we consider a holonomic property of such equations, valid in the most general 
situation, even without the assumption (2.5). 

3 Quasi-Holonomic Modules 

The class J-'n+i of functions, among which we looked for a solution of (2.3), is basic in the 
function field counterpart [15| of the theory of holonomic modules Let us recall its principal 
notions. 

Every operator a G 2t„+i can be written in a unique way as a finite sum 

A = Y,<^i,,,n,...,yd^K ■ • • An", ai,,,,„...,^ e K,. (3.1) 

We introduce a filtration in 2l„+i denoting by V^, z/ G Z+, the i^c- vector space of operators 
(3.1) with max{/ + ^ + ii + ■ ■ ■ + in} < v where the maximum is taken over all the non-zero 
terms contained in (3.1). Then 2t„+i becomes a filtered ring. 

Suppose M is a left module over 2t„+i with a filtration that is 

OJlo C OKi C . . . C M, M = y OJlj, 

i>o 

each 2Jtj is a finite-dimensional vector space over Kc, and Fj^OJlj C 97ljy+j for any v,j G Z_|_. 
The filtration is called good if the corresponding graded module is finitely generated. 

For a good filtration there exist a polynomial x ^ and a number G N, such that 

s 

dimmt, = J]dim(9Jli/9Jli_i) = x{s) ioi s > N 

i=0 

(dim means the dimension over Kc)- The number d{M) = degx, called the (Gelfand-Kirillov 
or Bernstein) dimension of M, does not depend on the choice of a good filtration. In general, 
d{M) < n + 2. In some cases d{M) can be evaluated explicitly. For example, d^f'n+i) = n+ 1, 
while (i(2t„+i) = n -|- 2, if is considered as a left module over itself. In contrast to 

the classical Bernstein theory (see [1]), in the positive characteristic situation d{M) can be 
arbitrarily small [15]. A module M is called quasi-holonomic if d{M) = n + 1 (in the survey 
papers [161 ITT] the term "holonomic" is used). As in the classical case, modules associated with 
basic special functions on K possess this property [15] . In general, the holonomic property is a 
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kind of a "quality certificate" to distinguish objects (equations, functions etc) with a reasonable 
behavior. 

If J is a non-zero left ideal in 2tn+i, then the left 2t„+i-module can be endowed with 

a good filtration, and 

d{%n+i/I)<n + l. (3.2) 
Returning to the equation (2.3) denote by R the operator in the left-hand side: 

i? = P(Ai,...,A„) + Q(Ai,...,A„)ci 

where P, Q are non-zero polynomials. Let I = St^+iP. 

Theorem 2. The module M = 2t„+i// is quasi-holonomic. 

Proof. Due to (3.2), we have to show only that d{M) >n + l. First we prove two lemmas. 
Lemma 1. An operator (3.1) is linear if and only if ai^^^i^^,,,^i^ = for I ^ fi. 

Proof. Let cr G Kc. Suppose that Aa = a A, that is 

$^o^a,,,,,...,„r'd^A- . . . Atr = Y.ai,,,,,...,y-\'d^A\^ . . . A^". 

By the uniqueness of the representation (3.1) [15], we find that cr^' = cr, whenever ai^^^i^^,,,^i^ ^ 
0. Since a is arbitrary, that is possible if and only if / = /i. ■ 

Lemma 2. The ideal I does not contain non-zero linear operators. 

Proof. Using the commutation relations (2.2) we can rewrite the representation (3.1) of an 
arbitrary operator A G in the form 

^ = E <M....n^ ■ ■ ■ ^nr'd^, <,,n,...,n e Kc (3.3) 

Just as in (3.1), the coefficients are determined in the unique way, and Lemma 1 

remains valid for the representation (3.3). 

Suppose that an operator A G 2t„+i is such that AR is linear. Let us write (3.3) in the form 

N N 
/=0 ^t=0 

where ai^ are the appropriate elements of the commutative i^c-algebra V (without zero divisors) 
generated by the linear operators Ai, . . . , A^. We have also R = j + Sd, ■jjS & V, ■y 0, 6 0. 
In this new notation, 

(N N \ 
5^5^a;^r'# (7 + W). 
1=0 ^=0 / 
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As an element 7 G is permuted with powers of r and d, in additional terms (appearing 
in accordance with (2.2)) the powers of r and d respectively are the same, while the degrees of 
elements from T> (as polynomials of Ai, . . . , A„) decrease. Therefore 

N N N N 

1=0 f_i=0 1=0 f_i=0 

where 7', 6' G 7' 7^ 0, 6' 7^ 0, whence 

N N N 
1=0 u=l 1=0 

By Lemma 1, 

aio = 0, / = !,..., AT; (3.4) 

aiN = 0, / = 0,l,...,iV. (3.5) 

Considering terms with I < N,!/ = N, we find that 

aiN'y' + a/,Ar-i5' = 0, 

and (3.5) yields ai^N-i = 0,0<1<N — 1. Repeating the reasoning we obtain that aiu = for 
l<u. 

On the other hand, for / > 2, i/ = 1 we get 

a/17' + CiioS' = 0, 

and, by (3.4), an = 0. Repeating we come to the conclusion that ai^, = for / > z/, so that 
A = 0. m 

Proof of Theorem 2 (continued). Let us consider the induced filtration in M (see [4J). The 
subspace DJlu is generated by images in M of the elements (3.1) with max(/+/i+2i + - ■ ■+in) < i^; 
those two elements whose difference belongs to I are identified. Let us consider elements with 
/ = /i. 

Elements of the form T^d^/S^^ ■ ■ ■ A^" with different collections of parameters . . ,in) 

are linearly independent in If some linear combination of their images equals zero in M, 

then the corresponding linear combination of the elements themselves must belong to /, which 
(by Lemma 2) is possible only if it is equal to zero. Therefore the images of the above elements 
are linearly independent, so that 

dimSOlj, > card{(/,ii,...,i„) G Z"+^ : 2/ + zi + ■ ■ ■ + 4 < z/} 

> card{(/,zi,...,z„) G Zl^^ : I + h + ■ ■ ■ + in < K2]} 

(in contrast to the rest of the paper, here [■] means the integer part of a real number). 

Evaluating the number of non-negative integral solutions of the above inequality as in the 
proof of Lemma 1 from [15j we find that 

dim 971, > f f''/^] + ^ + A > ci[z//2]"+i > C2Z/"+^ 
V n + 1 J 

for large values of u (ci,C2 are positive constants independent of z/). Thus, d{M) > n + 1, as 
desired. ■ 
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4 Hypergeometric Equations 



Let n > max(r, s) , r, s e N, 

r 



i=l 



Q{ti,...,tn)^ll{tj-bj), 

3=1 

where ai,bj e Kc, and the elements bj do not coincide with any of the elements [i/], u — 
0, 1, . . . , cxD. Then the condition (2.5) is satisfied, and the Cauchy problem for the equation 
(2.3) is wcll-poscd. 

Let us specify the initial condition (in terms of prescribing the values of Co,ii,...,i„) as follows: 

Co,0,...,0 = 1, C0,ii,...,i„ — 

for all other values of ii, . . . , i„. Then Cm,n,...,in — for all sets of indices (m, ii, . . . , i„) except 
those with m — ii — . . . — in- Denote — Cm,m,...,m.- By (2.6), we find that 



1=1 



(72 




and, by induction, after rearranging the factors we get 

fl([0]-a,)^'"([l]-a,)''"'-^--([m-l]-a,)^ 

Om^'-^ , m=l,2,.... (4.1) 

n([0]-y^-([l]-6,)^'"-^---([m-l]-6,)^ 

For a e Kc, denote (a)o = 1, 

{a)rn = m - a)^"([l] - a)^'"" ... ([m - 1] - a)\ m > 1 (4.2) 

(of course, [0] = 0, but we maintain the symbol [0] to have an orderly notation). The Pochham- 
mer type symbol satisfies the recurrence 

(a)^+i = ([m] - ay{a)l. (4.3) 

If a = [—a], q; e Z, then 

([m] - ay = (x'^"' - x'^-'^y = {x^'^^" - x) = [m + ar'"'"', 



so that in this case the recurrence (4.3) coincides with the recursive relation [19] for the 
Pochhammer-Thakur symbols (1.2). Our normalization (a)o = 1 is different from (1.2) and 
resembles the classical one. 

In the above situation, it follows from (4.1) that the solution u{z; ti, . . . , tn) of the Cauchy 
problem (2.3)-(2.4) is given by the formula 

u{z; ti, . . . , tn) = rFsiai, . . . , a^; 6i, . . . , 6^; ^1 ■ ■ ■ tnz) (4.4) 
where rFg is the new hypergeometric function 

OO / > , \ g-m 

^F,(ai, . . . , a^; 6i, . . . , 6,; 2) = 2^ — . • (4.5) 

m=0 

The series (4.5) has a positive radius of convergence since bi, . . . ,bs do not coincide with any 
of the elements [z/], = 0, 1, . . . , oo (such parameters will be called admissible) . 

If ttj = [— aj], bj = [—/3j] (z = 1, . . . , r; j = 1, . . . , s), G Z, j3j G N, then the function 
(4.5) coincides with Thakur's hypergeometric function (1.3) up to a change of variable z ^ pz 
where p depends on all the parameters but does not depend on z. 

Let hm be the coefficients of the power series (4.4), that is 

, _ (Ql)m • • • Mm _ n 1 9 

/h \ /h \ n ' ^ U, i, /, . . . . 

\"l/m ' ' ' \"s/ m-'-^m 

Since = [m+l] = (x'?'" - x"'")' = ([m] - [-1])", we find that 

hm+1 _ \ {[m] - ai) ■ ■ ■ {[m] - ar) Y ^^^^ 



hi l([m]-6i)---([m]-6,)([m]-[-l]) 



h 

The identity (4.6) means that the ratio — — is the g-th power of a rational function of [m], 

which is a clear analog of the basic property of the classical hypergeometric function. Note 
that any rational function of [m] may appear in (4.6), except those for which (4.6) does not 
make sense. Thakur's hypergeometric function corresponds to the case of rational functions 
with zeroes and poles of the form [i/] , z/ G Z. 

It follows from (4.4) that AjW = Au where A = rd is the difference operator (1.1) in the 
variable z. Therefore satisfies the equation 

n (A - a,) - - ^i) j ^1 rFs = 0. 

In particular, for the Gauss-like hypergeometric function 2-F1 we have 

{(A - a) (A - 6) - (A - c)d}2Fi = 0. 

Substituting A = rd and using the commutation relations (2.2) we can rewrite this equation 
in the form 

{r(l - T)d^ - (c - ([l]^/« + a + 6)r) d - ab} 2F1 = 0, 
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which corresponds to (1.5). 

If 6 7^ [u], 1/ = 0, 1, . . . , oo, then, in particular, |6 + x| > > 0, whence 

\b - = \ {b + x) - x''" \ = \b + x\> n 

for large values of v. This means that 

>/xi>0, z/ = 0,l,2,..., 

so that 

Therefore, if \z\ is small enough, then the series (4.5) converges uniformly with respect to the 
parameters G Kc and bj G Kc\{[z/], u = 0,1, . . . , oo} on any compact set. Thus, the function 
rFg is a locally analytic function of its parameters. 



5 Contiguous Relations 

Among many identities for classical hypergeometric functions, a special role belongs to relations 
between contiguous functions, that is the hypergeometric functions 2-^1 whose parameters differ 
by ±1 [1]. For Thakur's hypergeometric functions, analogs of the contiguous relations were 
found in [191 1201 [21] . Here we present analogs of the contiguous relations for our more general 
situation. Of 15 possible relations, we give, just as Thakur did, only two, which is sufficient 
to demonstrate specific features of the function field case; other relations can be obtained in a 
similar way. Note that the specializations of our identities for the case of parameters like [—a], 
a G Z, are slightly different from those in [191 ISOl [21] , due to a different normalization of our 
Pochhammer-type symbols. 

Our first task is to find an appropriate counterpart of the shift by 1 for parameters from 
Kc- Denote 

Ti(a) = (a- [l])i^ aeKc. (5.1) 

If a = [—a], a G Z, then 

m-a]) = {[-a] - [1])'/" = x^-"-' - X = [-a - 1], 

so that the transformation Ti indeed extends the unit shift of integers. The inverse transfor- 
mation is given by 

T_i(a) =a^ + [1], aG^c. (5.2) 

Theorem 3. The following identities for the Pochhammer-type symbol and the Gauss-type 
hypergeometric function are valid for every m G N and any admissible parameters from Kc: 

(Ti(a))m = a~'''"(a - [m])(a)m, a 7^ 0; (5.3) 
(a)„+i = -a^™^'(Ti(a))^; (5.4) 
(T-i(a))„ = -([l]+aT'"(a)L_,; (5.5) 

{T.,{a))m = -7^^^f^(«)- (5-6) 
([m — IJ — aj'J 



11 



2Fi{Ti{a),b;c;az) - 2Fi{a,Ti{b); c;bz) = (a - 6)2^1(0, 6; c; 2); 



(5.7) 



2Fi(a, b; c; ;2) - 2^1(0, 6; c; ;2)« + (c^ - 6^)2^1 (a, 6; Ti(c); c-i^)« 

- (a« + [l])2Fi(T_i(a),6;c; (a'? + [1])-^^) = 0. (5.8) 

Proof. Substituting (5.1) into (4.2) and using the fact that H + [l]^/^ = x"" -x^^" = [z/+l]^/'', 
we get 

(Ti(a))„ = ([1] - ay"'-\[2] - a)^'""' ■ ■ ■ ([m] - a), (5.9) 

which imphes (5.3). If we raise both sides of (5.9) to the power q and compare the resulting 
identity with (4.3), we come to (5.4). The proofs of (5.5) and (5.6) are similar, based on the 
identity [u] - [1] = [i/ - 1]<?. 
Using (5.3) we find that 

2Fi(Ti(a), b; c; az) = ^(a - [m])^^'^'^^^'^"' 



(c)m-D. 

2F,{a,T^{b)■c■bz) = ^(6- [m]) 



m=0 



m=0 

which implies (5.7). Similarly, if we write down all the terms involved in (5.8) and use the 
identities (4.3), (5.3), (5.4), and (5.6), after rather lengthy but quite elementary calculations 
we verify the required identity (5.8). ■ 
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